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Abstract 

We discuss the convergence of cosmological perturbation theory. We 
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jrt ' prove that the polynomial enhancement of the non-linear corrections 

expected from the effects of soft modes is absent in equal-time cor- 
relators like the power or bispectrum. We first show this at leading 
^ ■ order by resumming the most important corrections of soft modes to 

an arbitrary skeleton of hard fluctuations. We derive the same result 

l/~) I in the eikonal approximation, which also allows us to show the absence 

of enhancement at any order. We complement the proof by an explicit 

■^ , calculation of the power spectrum at two-loop order, and by further 

numerical checks at higher orders. Using these insights, we argue that 
the modification of the power spectrum from soft modes corresponds 
at most to logarithmic corrections. Finally, we discuss the asymptotic 
behavior in the large and small momentum regimes and identify the 
expansion parameter pertinent to non-linear corrections. 



1 Introduction 

The theory of cosmological perturbations is the standard tool to understand 
the emergence of the large scale structure of the Universe [H [2] . This ap- 
proach is based on the assumption that small perturbations around an other- 
wise homogeneous and isotropic Universe grow with time due to their gravi- 
tational interaction. This growth is particularly efficient for scales inside the 
Hubble horizon and in the matter-dominated epoch [1]. The tiny amplitude 
of the primordial perturbations allows for a perturbative treatment of this 
non-linear problem, but the aforementioned growth eventually invalidates 
this approach. The scale at which this happens has been traditionally iden- 
tified with moments of the linear power spectrum. A particularly interesting 
scale is created by the enhancement of non-linearity coming from the cou- 
pling of modes of small momenta (soft) to the modes of the scale of interest 
{hard), 

k-^\^ [dqP'^iq,^), (1) 



where P^{q,ri) is the linear power spectrum at time t] (we will be more 
precise about this point below). To study physics beyond this scale, it seems 
unavoidable to resort to non-perturbative schemes. Particularly well-suited 
to deal with this effect is the scheme known as renormalized cosmological 
perturbation theory (RPT) P]. This systematic approach to cosmological 
perturbation theory is based on the introduction of the non-linear propagator, 
for which the non-linear effects associated to the scale ([1]) were derived in [1]. 
These result in an exponential suppression controlled by kNL that was also 
conjectured to be the leading effect in the matter power spectrum. 

The seminal work [3] was followed by many other resummation schemes 
put forward to cope with non-linear effects in an analytical way (i.e. without 
resorting to large A^-body simulations) [5l |6l [TJ El [9l [10] . The hope is that 
after resummation of a subset of diagrams the resulting expansion is under 
better control. For schemes related to the effects of the scale ([1]) in the power 
spectrum this expectation is at odds with two results. First, it is known since 
a long time that the leading contributions to the power spectrum from soft 
modes at arbitrary loop order cancel fiT\ [T2] . In the above mentioned re- 
summation schemes this cancellation is not explicit. Second, one can try to 
systematically understand the effects of soft modes on hard modes by using 
the eikonal approximation [13]. In this case, it was shown that the men- 
tioned suppression for the propagator is present while the power spectrum is 



unchanged [H] . A first aim of this work is to reconcile both approaches and 
understand how the eikonal result can be recovered from the diagrammatic 
technique of resummation of soft modes. 

The cancellation of the effects from soft modes suggests that the conver- 
gence properties of standard perturbation theory (SPT) are not ameliorated 
by resuming the soft modes. It also suggests that the effects in the power 
spectrum associated to k^L are spurious, which immediately rises the ques- 
tion of what is the real parameter governing the non-linear corrections for 
this observable, also in SPT. It is simpler to answer this point after imple- 
menting the eikonal approximation in a controllable way. We will do this 
in the second part of the paper. The scale (or parameter) controlling the 
non-linear dynamics in this case, and not ([T]), should be the one tamed by 
any resummation scheme with a better validity than the standard case for 
the non-linear dynamics. 

Our work is organized as follows: Section |2] clarifies our notation and 
reviews standard perturbation theory (SPT). Next, section E] discusses the 
two main resummation schemes we are concerned with, namely renormalized 
perturbation theory (RPT) and the eikonal approximation. In section H] we 
generalize the resummation of RPT to a larger class of soft corrections that 
are relevant to the power spectrum and other equal-time correlators. This 
motivates the discussion of next-to-leading order corrections in the eikonal 
approximation presented in section El The main result is that no enhance- 
ment from soft vertices should be present in equal-time correlators. Finally, 
we further support this claims by explicit analytic and numerical results pre- 
sented in section |6] before we conclude in section [71 Some technical details 
are relegated to the Appendices. 

2 Standard cosmological perturbation theory 



In this section we set up the system of equations relevant for our discussion. 
As stated in the introduction, we are interested in understanding some fea- 
tures about the behavior of cosmological perturbations when their amplitude 
grows to the point that non-linear corrections are important. For this prob- 
lem it is enough to consider sub-horizon perturbations in a matter-dominated 



era. We will also assumqj that the matter in the Universe is a perfect fluid 
described by a density field p{x,t) = p(t)(l + 6{x,t)) and a velocity field 
v^{x,t) defined at a conformal time r = J dt/a{t). In this case the New- 
tonian cosmological perturbations yield an accurate description [1]. Finally, 
we will reduce our analysis to the case without vorticity where 6 = diV^. In 
this case, the perturbations can be written in a two-component form [21 [T7] . 

d^n^aik, 7]) + ilab^f^ik, T]) = -fabcih, k2)^b{ki, v)^c{k2, v) , (2) 

where 

<iJ,ik,v)^Sik,v), <iJ2ik,v)^-^^, (3) 

and we have introduced the functions V. = ^J^'^ and f+{r]) = ^"na( ) ^^"^ 
the time i] = lnD_|_(r), with D^{t) being the linear growing mode of the 
density contrast 

S'^{k,T) = D+ir)6oik). (4) 

The momenta k, ki, ^2 are vectors and the equation is understood as being 
summed over double indices and integrated over the two momenta ki and 
k2 with the measure 6^^\k — ki — k2) on the right-hand side. This fixes our 
convention for Fourier transforms. In the following we only write indices and 
integrations when the notation of an equation could be ambiguous. 

The matrix 7 constituting the mode coupling can be written in symmetric 
form with the elements 

'li2i=a{ki,k2)/2, -fn2 = a{k2,ki)/2, 7222 = /3(A;i, fe), (5) 

with 



{ki + k2) ■ ki 

If 
{ki + k2fki ■ k2 

01,2 i,2 



a{k,M = ^J^LZ^IL^^ (6) 

f3{kuk2) = ^^-2T2 , (7) 



and all other elements vanishing. 



^Deviations from this assumption can be taken into account by using the effective 
language of [15] (see also [IS]). Since we focus on the first non-linear effects those deviations 
are not important for our results. 



In the case where one of the momenta flowing into the vertex is soft 
{ki <^ k OT k2 "^ k) the previous vertex reduces to 

lijk -^ lijk = Oj2dik^77a- + ^fcsdii^jp- • (8) 

The matrix Q depends on the underlying cosmology and can depend on the 
conformal time rj (but not on momentum). It may also contain information 
about modifications of gravity [9]. Remarkably, our main results will be valid 
for any Q. For completeness let us remind that in a fiat matter-dominated 
Einstein-de Sitter Universe for which Q^ = 1 and D^{t) = a{T) it reads 

"=(-3/2172)- w 

Note also that the common Zel'dovich approximation [18] differs from the 
exact dynamics ([2]) only by the choice of the r2-matrix [21 [191 120] 



^"^ = (0 -\) ■ ^^°) 

The equations (Ej) can be solved perturbatively by treating the right side of 
the equation that mixes modes with different momentum as a perturbation 
to the linear equation 

9^*^(A;,r/) + fi(r/)*^(A;,r/) = 0. (11) 

This perturbative scheme, known as standard perturbation theory or SPT, is 
based on the assumption that the amplitude of the perturbations \E' is small. 
The solution to the previous equation can be easily written in term of the 
Green's function and the initial conditions \E'(/c,?7o) as 

^\k, v) = e- ^''0 '^'^^(^ ^(fc, r]o)e{v - Vo) = 9{v, ^o)*(fc, Vo) , (12) 

where we used the notation 0(x) for the Heaviside step function to avoid 
confusion with the velocity field 6. If the matrix ^2 is independent of rj, the 
Green's function depends only on the difference rj — tjq. For the case iQ, its 
explicit form is 

g(»)-w) /3 2\ e-3{»7-»7o)/2 / 2 _2\ 
9{v,Vo) = — ^ — (3 2)®(^~^o) + ^ (3 3 ) 0(^-^0), (13) 



from where we can readily identify the growing and the decaying mode. The 
solution to the equation (j2]) can formally be written 

^aik,r]) = gabiv^Vo)^bik,rio) 
rv 
+ 9ab{r].n)lbcd{Kk2)^c{Kfi)^d{k2,f])- (14) 

Jrio 

The perturbative solution (in powers of \l/a(/c,?7o)) can then be obtained by 
successively reinserting the left-hand side into the mode function 

^a(fc,^) = gabiv,Vo)'^bik,r]o) 
rv 
+ / df]g{r],f])ab1bcdiki,k2) [^ce(r/,r/o)^e(A;i,r/o)] [^d/(r/,r/o)*/(A;2,^o)] 

JriO 

+ ■■■. (15) 

Physical observables can be evaluated by performing statistical averages over 
these fields. A particularly useful one is the power spectrum given by 

{^a{ki,v)^b{k2,r]))=6^'\k, + k2)Pab{kuv)- (16) 

The conservation of momentum in the mode coupling term ensures that the 
two point correlator (II 6p is proportional to the delta function at all times. 
This reflects the translation invariance of the system. The leading order 
contribution to it in SPT is the so-called linear power spectrum 

Pabiv) = 9aciv,Vo)9bdiv,Vo)Pcdik). (17) 

where we have introduced the initial power spectrum 

{^aik,,Vo)'^,{k2,Vo))^5^'\k, + k2)P!^,{k,). (18) 

Thus, 

Pabik,rj) = P^^{rj) + .... (19) 

In the following we only deal with Gaussian initial conditions that are uniquely 
specified by the initial power spectrum P^{k). 

The perturbative expansion can be pictorially represented by classical 
Feynman rules [3] (see also [211 [22] for the similar treatment of the fluid 
equations in the context of turbulence). The two building blocks are the 
vertex and the Green's function displayed in Fig. [1] The causal structure of 
the integral flT5|l enforces a flow of time that is indicated by the arrow. The 
initial power spectrum is indicated by a box. The leading topologies that 
contribute to the power spectrum are shown in Fig. [21 
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9ab{V, V) y ^ , . Pah{k) 

(T],a) {fi,b) (/i;i + 4,a)\^"'^ ' (k, a) i-k,b) 

{k2,c) 

Figure 1: Building blocks of the Feynman rules of standard perturbation theory. 
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Figure 2: All tree level and one-loop contributions to the power spectrum. 

3 Renormalized perturbation theory and the 
eikonal approximation 

The expansion in loops is reasonable even though the equations ([2]) do not 
have an intrinsic small expansion parameter, a coupling constant for instance. 
This is because every loop comes with an additional initial power spectrum 
in SPT. This can serve as an expansion parameter as long as the power 
spectrum is small in adequate units. 

Even in the limit of a small initial power spectrum, the convergence of 
the perturbation theory is not guaranteed. The first reason is obvious from 
Eq. (TT5!) : the initial power spectrum appears in the perturbative expansion 
multiplied by the Green's function. Since the latter contains a growing mode, 
perturbation theory tends to be worse behave as time grows. Another effect 
that can potentially lead to the breakdown of perturbation theory is related 
to the enhancement of the vertex (coupling) between hard and soft modes. 
For example, the diagram in Fig. [3] scales for k ^ q a.s 

r r ~ ~ f {k ■ qY 

9abiv,Vo) df] dr]g2cif],Vo)92div^Vo) d^q ^—Pcdil)- (20) 

So even if the final integral is finite and the growing of the Green's function 
small, there results an enhancement from attaching more soft loops to a line 




(/c, a) k + a (k,b) 



Figure 3: One- loop contribution to the propagator, 
with momentum /c if /c is large enough, i.e. ii kad^ 1, where 

aj{A,v) ^^f£q^-^PHQ,v) = ^ f dqP\q,r^) , (21) 

where P^ = P22 and A represents the scale at which we cut-off the initial 
power spectrum. This is precisely the scale of non-linearity discussed in the 
introduction, Eq. ([1]). In particular, /c^cr^ can be larger than the dimensionless 
variance of the density field [12] 



.A 

afiA,r])=A7Tj dqq^P'^^q.r,), (22) 

which would control the non-linear scale in the absence of vertex enhance- 
mencl. At this stage it seems that to be able to describe physics beyond the 
scale ad one should first be able to sum up all soft contributions to the prop- 
agation of hard modes. Once this is done, one expects to find a perturbation 
theory with better convergence properties. This idea was put in a firm basis 
in the theory of renormalized cosmological perturbations (RPT) [31 H] (see 
also [211 122] for similar resummations in the case of fluid dynamics). In this 
case, the role of the Green's function (TT3l) in the Feynman rules is played by 
the propagator 

Gi;.(M,.).-(.-.o.(^f;^). (23) 

that in terms of perturbation theory contains all diagrams with exactly one 
incoming and one outgoing line and arbitrary (soft or not) loop corrections 
to it |. 



^We chose the notation for these two quantities that best adapts to previous literature. 
We also define crdj(A) = (Td,i{A,ri = 0) for the values today, and ad = crd{A — > 00); 
note that this limit exists for realistic power spectra. In contrast, (Ji{A) has a logarithmic 
sensitivity to A, as will be discussed in detail in section [71 

•^The meaning of the superscript (1) is clarified below. 



As we will review in the next section, the leading behavior for large mo- 
menta k (and considering only the contribution from the growing mode) can 
be resummed for the propagator of renormalized perturbation theory [3] 

GRPTik,V,v) =^ab(^,^)exp i--k^al{a{r]) -a{f])y\ , (24) 

where we used 0(77) ^ e^ and set a = 1 today. Even if the vertex is not 
modified in this framework, the steep suppression of the propagator at large k, 
in principle compensates the enhancement provided by soft loops attached to 
this propagator, which makes one expect a better convergence of perturbation 
theory in this case, with the non-linear scale set beyond ad- 

The perturbative expansion in RPT can be formulated in terms of the 
n-point propagators defined as 



Gi:!...„ (k, ,, ,mk - E *.) - ^ ( ... J:!:^'^. a_ r.) ■ (2«) 



In terms of those, the power spectrum at late times can be written as 

n \i=l ^ J 

X G^:^l„,Jk,v,Vo)Glt.b^-k^,V,Vo) . (26) 

Pictorially, C^") is the sum of all diagrams with n incoming lines and one 
outgoing line. As long as the resummed n-point functions G^"^ are well 
behaved for large k, the sum should converge. It was shown in [23] that 
the corrections from soft modes to the individual n-point propagators also 
generate an exponential suppression related to the scale cr^. 

For the power spectrum and related observables, the previous results 
provide only a partial resummation of the effects of soft modes. This is 
commonly encoded in an expression of the form (compare with (fT9l) ) 

P = {c^nl^Ty P' + Pmc, (27) 

where the piece Pmc contains all the contributions from other diagrams. 
From the previous arguments, one may expect Pmc to be better behaved 
than the corrections in Eq. (TT9|) and that RPT with the leading effect of the 
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soft modes to the propagator resummed have better convergence properties 
than SPT. 

To check this, one needs to deal with all soft corrections to the hard 
skeletons behind observables. This is the purpose of the eikonal approxima- 
tion [13]. The main observation of this approach is that in ([2]) the contribu- 
tion from the soft modes to the mode coupling is approximated by 



d^hd^k2 -fabcih, h)'^b{ki, r])<i'^{k2,v)S^^\k -h- k^) 



Using this relation in (|2]) leads to the result (cf. f lT2|) ) 



'^{k,7]) = g{r],r]o)exp 



dt] I d^q^—'$2{q,v) 



Vo 



^(A;,7/o) 



(28) 
(29) 

(30) 



In turn, for the propagator (123|) in the eikonal approximation this yields 
immediately 






eikonal 



{k,V,v) = ^(^,^)(exp 



73..^-?, 



dfj / rf^g— ^^2(g,^) 



q- 



(31) 



that for Gaussian initial conditions is (in leading order, see (fT2|) ) given by 
the second cumulant 



G'ilonJk, V, V) = giv, V) exp (^) 



(32) 



with 



C2 



dt] / d^q^^'$2{q,v) 

fj J H 

V 



^- —^2yq,''l)j I 



connected 



dfidrig2a{fl,Vo)92b{'n,Vo) d q 



ik-q)\o 



(33) 



PM 



To compare it with (124|) . we consider only the growing mode, for which 

C2 = -k'al{a{r^)-a{r^))\ (34) 

which shows that the eikonal approximation indeed resums the leading soft 
corrections. A possible physical interpretation of this propagator is that the 
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Figure 4: Two-loop contributions to the propagator. 

modes of hard momentum k are scattered by the background of soft modes. 
This decorrelates the modes over time what leads to the exponential fall-off 
of the propagator for late times. 

One of the advantages of the eikonal approximation is that it allows to go 
beyond the propagator and compute the effects of the soft modes in the power 
spectrum [TH]. Since the exponent in (pi]) is odd under a sign flip of k and the 
power spectrum involves the combination (\i/(/c)\E^ (/;;')) oc S{k + k'), the two 
exponential factors cancel and the resulting power spectrum is unsuppressed. 
Hence, the complete resummation of the leading effect of soft modes does not 
produce an expression as fl27|) (see also (126|) ). but leaves flT9|) unchanged [T3] . 
This result can also be derived in the diagrammatic language of RPT, which 
we do in the next section. 



4 Resummation of general skeleton diagrams 

Let us start this section by reviewing the resummation of the leading contri- 
butions to the propagator as first given in [3] . We want to resum all leading 
soft corrections in the limit of large k. A vertex with one soft and one hard 
incoming mode is enhanced by a factor k ■ q/q"^, while a vertex with two 
similar modes is not enhanced. Hence, for fixed loop order the dominant 
contribution comes from attaching all soft modes directly to the hard mode 
that flows through the linearized propagator. Possible diagrams at one and 
two loop order are shown in Figs. |3] and HI In the limit of soft corrections, 
the vertex is proportional to the Kronecker delta ([8]). The inflowing momen- 
tum can be neglected at leading order and the Green's functions 5'(?7i+i,''7i) 
along the hard mode combine to one linearized propagator involving only 
the very first and the very last time, 5'(?7,''7o)- This is a generic result that 
hinges on the group structure of the propagators. In particular it is valid 
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k,V 



k 



O 

01 Qo Q2n-lCl0 



k,Vo 



Figure 5: Soft modes attached to a hard hnearized propagator. If the inflow of 
momentum is neglected, the order of the vertices is irrelevant. 

for any expansion history of the Universe and beyond the restriction to the 
growing mode case. A diagram with n loops involves attaching 2n linearized 
propagators of the form depicted in Fig. O As long as the inflow of soft 
momentum is neglected, the order of all the vertices is irrelevant. The time 
integrations involves (with rj2n+i = v) 



2n 



Vi+i 



Yl I driig2a{vuVo) 



(35) 



Since we are considering the effect of all possible contractions with the initial 
power spectrum, the symmetries of the resulting expressions allow us to 
rewrite the time integrations as 



1 ^"' f^ 

^n / driig2a{Vi,Vo)- 



(2n)! 



(36) 



i=0 -^ 'JO 



In total, there are (2n— 1)!! possibilities to contract the soft modes with initial 
power spectra. Thus, at n-loop order one finds a contribution (c2/2)"'/n! and 
summing over loops the eikonal result fj32|) is recovered. Notice that this 
result is valid for any matrix Q, which just changes the form of the linearized 
propagator and hence the particular value of ad- 

4.1 Power Spectrum 

Next, we consider all soft corrections to the power spectrum. A generic graph 
at four-loop order is shown in Fig. [61 Imagine there are riu loop corrections 
to the left hard linearized propagator and rirr loop corrections to the right 
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Figure 6: A contribution to the power spectrum at four-loop order with nu = 2, 
nir = 1 and rirr = 1- 

linearized hard propagator. In addition, we denote the number of soft loop 
corrections connecting the left linearized hard propagator with the right one 
by nir. Following the same logic that applied for the propagator, one can 
extend all time integrations from rjo to t] what leads to a factor 1 / {2nii + nir)l 
for the left linearized hard propagator and l/{2nrr + nir)\ for the right one. 
There are ( ^"^"^'r^ combinations to split the soft modes connected to the 
left branch into the two groups, such that the combinatorial factor before 
contraction is l/{2nii)\nirl and similarly for the right soft modes one finds 
1 / {2nrr)lnir\ ■ As before, there are (2nii — 1)!! possibilities to contract the 
left loops and (2n,.r — 1)!! to contract the right ones. In addition there are 
riirl combinations to contract the left modes with the right ones. Finally, 
notice that the left and right loops produce a factor C2 while the soft loops 
connecting the left branch with the right one leads to a factor — C2. In 
summary, one finds that the leading corrections from soft modes for the 
power spectrum of a hard mode reduce to the factor 

nii,nir,nrr 

Thus, the corrections that connect the different branches exactly compensates 
for the exponential suppression found in resumming the soft corrections to 
the propagators. This is in complete accord with the eikonal approximation 
that predicts a suppression in the propagator but not in the power spectrum. 
This cancellation of the leading soft corrections to the power spectrum was 
also observed in [TT] . 

The previous result implies that the expansion fl26l) does not necessarily 
improve the convergence in the power spectrum compared to standard per- 
turbation theory. According to the above resummation, in the limit where all 
ki but one momentum are soft the ra-point functions C-"^ are enhanced such 
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that the sum produces an exponential that cancels the exponential suppres- 
sion in each of the G^"-* observed in RPT [23]. Parametrically, the leading 
soft corrections to the different contributions to the sum fl26|) with resummed 
n-point propagators scale as P ~ Sn ^(^'^d('7))^'^exp(— /c^o"d(?7)^). So for 
large /c^0"d(?7)^ ^ 1, the sum can only start converging after n ~ k^crdirj)"^ 
terms are taken into account. In fact the situation is even worse, since re- 
producing the leading /c-behavior order by order is not enough to ensure the 
cancellation of subleading terms in the regime of large k as we show below. 

4.2 General Skeleton 

The soft corrections to arbitrary skeletons of hard modes can also be re- 
summed. For example consider any of the one-loop contributions to the 
power spectrum depicted in Fig. |2] in a regime where none of the involved 
momenta is soft. The diagram involves several linearized propagators with 
momenta fcj. Soft loops can connect arbitrary linearized propagators and 
we denote the according number as Uij. The total number of soft lines at 
a linearized propagator i is then Ni = 2njj -|- J2j^i''^ij- -^^ before, the soft 
lines are time ordered and extending the time integrations over the full range 
leads to a factor 1/Ni\. The full range is hereby given by the time the original 
linearized propagator in the skeleton diagram depends on that we denote rji 
and f]i. Splitting the soft lines into groups yields several binomial factors 
that contribute factors of the form l/(2njj)!/ ]^ . /^rzjj!. This is valid for an 
arbitrary skeleton of hard modes. 

Connecting a soft line of a linearized propagator i with a soft line of a 
linearized propagator j gives the integral 

- /rf^g^^P^Vg), (38) 



and the time dependence! 



q- q- 



Vi 



dv92a{v,Vo) / d7ig2b{fi,r]o). (39) 



^ Note that, since time has to increase along the arrows of each Unearized propagator, 
cf. (|13p , the time 77 of a soft vertex that is attached between vertices at fji and r]i has to 
He within the time interval fji < f/ < rji. 
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k2,Vo 



ki,Vo 



Figure 7: An example of several contributions to Yab that can be combined. 

With the corresponding combinatorial factor to connect soft lines, the final 
correction to the diagram reads 



exp 



-- / d'qY,,Y,,PM 



with 



Y„, 



E 



Vi _ ki ■ a 
dfj-^gab{fi,rio) . 

i. ^ 



(40) 



(41) 



This result is actually expected from the eikonal approximation as we will 
see in the next section. 

Depending on the skeleton diagram, the function Yab can be further sim- 
plified. Consider two linearized propagators that originate from two initial 
power spectra and end in a common vertex. The two corresponding con- 
tributions to Yab then share the same range of time integration and can be 
combined to one contribution that is proportional to the sum of the two mo- 
menta of the skeleton diagram. Hence, this term can also be combined with 
the subsequent propagator. An illustration of this is given in Fig. [71 This 
fragment of a diagram acquires through soft corrections the contribution 

Yab 3 / dfj — —gab{v,Vo)+ dfj — —gab{v,Vo) 

Jr]o q Jrjo 1 

+ / dri — gab{v,Vo) 

Jrji 



"'Ah + htq ,. . 

dr] ^ gabiV, Vo)- 

vo ^ 



(42) 



Remarkably, the functions Yab vanish for any equal-time correlator as for 
example the power spectrum and the bispectrum. 
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We will see in the next section that this will allow us to go beyond the 
leading order considered in previous analysis [HI [T3] and show that the effects 
of soft modes over hard modes k produce at most an enhancement of log k 
compared to the linear power spectrum. 

5 Eikonal approximation - the higher orders 

In this section we extend the previous results and show that also subleading 
soft corrections cancel each other in equal-time correlators. This is done 
using the language of the eikonal approximation [TTl [T3] . 

We want to separate in a systematic way the contribution in ([2]) that 
incorporates the enhancement from soft modes and that is used to resum 
the leading term in the previous analysis. To this aim we add and subtract 
several terms to the right-hand side of ([2]), and rewrite it in the equivalent 
forno 

2F,(/ci, k)%hciki, /cs) ^fe(/ci, r/)^,(/c, r/) 

+2F,{ku k)%i,ciki, ki) ^b{ki,r,){^,{k2, v) - ^dk, r/)) 

F^ki, k) {-fabc{ki, k2) - %bc{ki, ^2)) ^b(A;i, ?7)^c(A;2, v) 

Fe{k2, k) {'Jabciki, /cs) - %cb{k2, ki)) ^b(A;i, r])'^c{k2, v) 

+ (1 - F,(A;i, k) - F,{k2, k))-iabc{ku k2)^b{ku v)^c{k2. v) , (43) 

where we use the notation 

ki ■ (fci + ^2) 

The function F^ is a filter that distinguishes soft modes from hard modes, 
for example the Heaviside step-function F^ki.k) = Q{e\k\ — \ki\). In the 
previous construction only the first term on the right side becomes large in 
the limit \ki\ <^ \k\. This term can be resummed and absorbed into the 
Green's function 



7abc{kl, k2) = 5b25ac ^ W2 ' (^4) 



geikonai{k,'n,r]o) = g{7],r]o)exp 



dfj d'^qF,{q,k)—^^2iq,v) 



9{v,Vo)^ik,r],r]Q). (45) 



^Recall that a term 6^^'{ki + k2 — k) and integrations in fci and k2 are present in all 
the non-linear terms. 
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The four remaining terms in ( 143 p can be treated perturbatively. The 
calculation of averaged quantities at a given order involves mixed cumulants 
of the field due to the exponentiation of soft fields in fHSjl . At leading order 
in the soft background field, these cumulants reproduce the Y functions ( HTj) 
appearing in ( HOj) [T3] . 

The perturbative expansion defined out of ( 1431) with the propagator f HS]) 
seems better behaved than the one of SPT. This is because it tames the 
effects related to the scale (12T|) . This advantage would justify its use despite 
its higher degree of complexity coming from the proliferation of verticeqj 
and the appearance of mixed cumulants mentioned above. But the results 
from the previous section show that the leading effect in a^ cancels for equal- 
time correlators and suggest that this may be a spurious scale. We are now 
going to prove this for all subleading effects. Thus, for these observables the 
convergence of any scheme based on resumming the effects related to fj2T|) 
cannot be better than SPT. 

We support this claim by interpreting the split of fH3|) in terms of di- 
agrams. In any diagram, the putative enhancement at scales ka^ ^ 1 is 
related to soft momenta dressing a hard skeleton with momenta of order k. 
In the previous section, we showed that the leading effect (which at n-loop 
order scale as (fccTd)^"') cancels. To prove the cancellation of the subleading 
terms down to order 0{k^), we need to use the perturbation theory based 
on (143 p . In comparison to the eikonal limit this requires to account for the 
momentum injected into the diagram from the soft loops, for the full vertex 
compared to dH]) and finally for the self-coupling of soft modes (i.e. soft 
modes coupled to other soft modes). Notice first that the vertices coming 
from the last three terms in fH3|) do not produce any enhancement and can 
be treated perturbatively. The self-coupling of the soft modes is in fll^ rep- 
resented by the fact that we resum the time-dependent contribution \l/(g,?7) 
(the whole non-linear form). Whenever this propagator is used for vertices 
that conserve momentum, the phase will disappear. Thus, the only place 
where the scale ad may appear is in the term that does not conserve momen- 
tum, 

2F,{k^, k)%Uki, k2) ^b{ki, r/)(*,(fc2, r/) - ^,{k, 7])). (46) 



^ In particular some of the terms that have been added and subtracted in (|43l) are not 
translation invariant, which leads to an apparent violation of momentum conservation in 
individual vertices and a perturbative expansion where single contributions to the power 
spectrum are non-diagonal in momentum space (although the full theory still respects 
translation invariance). 
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li ki <^ k, the vertex is enhanced but the contribution conies with a factoiu 

(^(^2, V) - '^{k, 7])) ~ g{ri, ?7o)^(/c2, V, Vo) 

X [^{k2, Vo) - ^(A;, VoMk - /C2, V, Vo)] ■ (47) 

For equal-time correlators the factor C,{k2, tj, rjo) will combine with the remain- 
ing factors of the diagram to yield unity, while the second term in brackets 
is not enhanced in k any more and can be treated perturbatively (recall the 
factor F,(A;i, A;) in (USD) 

e(/c - /C2, V, m) = ^{h, r],Vo)-l+ / d'q ^^2(9, v) ■ (48) 

Hence, {'^{k2,f]) — "^{k,!])) is of order ki/k, which cancels the enhancement 
coming from 7, and the factors ^{k, r], rjo) also cancel in this vertex. 

This proofs that the enhancement from soft modes characterized by the 
scale ad is completely absent in equal-time observables. This includes the 
matter power spectrum and any other n-point correlation function, like the 
bispectrum. Notice that the argument in this section is rather general. 
The main assumption is that the leading term in f H3|l can be treated non- 
perturbatively and that the resulting cumulants are not enhanced for large 
k. Another assumption is that the Gaussian random field "^{k) can be ex- 
panded in a Taylor series, see also [H]. On the other hand, the precise form 
of the vertices is irrelevant for the argument to work as far as the full vertices 
respect translation invariance and provided that the soft-enhanced terms are 
proportional to the unit matrix '^abc ot 6ac such that ^ has no matrix struc- 
ture. The latter condition is for example violated in a multi-component fluid 
for the non-adiabatic decaying isodensity modes [11]. Notice also that, since 
the proof does not rely on the form of the linear part Q (given by (l9l) for Ein- 
stein de-Sitter), it is true for any cosmological expansion history, including 
any possible term that modifies the linear term. This means that it is valid 
as well for the Zel'dovich approximation f lTO|) . in accordance with [T2] . 

In the next section we complement the general proof presented above by 
an explicit check of the cancellation for the matter power spectrum in an 
Einstein-de Sitter cosmology at two-loop order, and for the loop integrand 
up to four-loop order. 



''This expression is precise to order 0{e). 



6 Matter power spectrum 

In this section we want to clarify the consequences of the cancellation of 
subleading corrections related to the scale ad for the matter power spectrum. 
The absence of physical effects related to the scale cr^ was already derived at 
leading order and all loops for the matter power spectrum in [H]. This was 
interpreted in [12] as a consequence of Galilean invariance. In [T21, [23] the 
explicit calculation for the matter power spectrum and bispectrum to two 
loops and with scale invariant initial power spectrum was performed, and no 
enhancement was found (the corrections are related to the scale cx;). Note 
that for scale-invariant initial spectra with a power-law index smaller than 
—1 the cancellation of soft corrections is closely related to the cancellation 
of infrared divergences [12] . 

Even if our proof of section [5] for the absence of enhancement by soft 
modes is valid for any initial power spectrum, in this section we will work 
with a realistic ACDM form. For this, an analytic expression is given by the 
prominent fit by Eisenstein and Hu [2S]. For analyzing the leading asymptotic 
behavior of loop corrections, it is often sufficient to consider the restriction 
of this fit to the Einstein-de Sitter case with spectral index Ug = 1, 

with 
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L^He + lMm, C^^^-^+TTeoim^'' /3-6.86MPC. (50) 

For the normalization of the spectrum we fit the previous formula to the 
linear power spectrum at 2; = found in [26], which yields a ^ 17000 Mpc^. 
Note that, for our numerical analysis, we use as input a linear power spectrum 
obtained with CAMB for WMAP5 parameters as provided in [26] . 

6.1 Formalism 

The different contributions to the power spectrum in standard perturbation 
theory at n-loop order can be written schematically as 

1+1 „ 
Pn~ioopik, V) = ^ d^h . . . d^kn+i A^^/ci, . . . , kn+i) 

m=l -^ 

xP\k^,r,) . . . P'^ikn+i, V) S^'\k -k,----kj. (51) 
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For example at one-loop level, the two diagrams in Fig. |2] involve the spectral 
dependence P^{k) P^{q) and P^{k — q) P^{q), which corresponds to m = 1 
and m = 2, respectively (as well as g = k2)- Those diagrams are related to 
the terms denoted bja P13 and P22 in the standard perturbation theory [2]. 
The expression of Am in terms of the symmetrized kernels F„(/ci, . . . , fc„) 
that characterize the non-linear evolution of the density field in the growing 
mode in an Einstein-de Sitter Universe [H [211 EZ] is shown in Eq. f l80p . 

We are interested in the case where the momentum k is much larger than 
the scale aj^. In this case, the previous integral seems to be dominated by the 
regime in which all but one of the modes ki are soft. This happens because the 
Am are homogeneous rational functions in the momenta and in the previous 
regime the largest enhancement is of order Am oc Ylii^ ' ^i/^iY "^ |A;p". 
Furthermore, the steep fall-off of the power spectrum fll^ at high k implies 
that the regime kaa^ <yi{k) exists. Thus, we may concentrate on the non- 
linearities associated to ad- The results in the previous section imply that 
this is a naive expectation: there must be a cancellation between different 
terms that eliminate this effect. To find this cancellation at different loop 
orders, we first use the Dirac-delta function to perform the integration over ki 
in (1511) . For large /c, enhanced contributions from soft modes would originate 
from the domain \ki\ <^ \k\ for i = 2, .., n -|- 1 and ki ~ k. In order to isolate 
these terms, we Taylor expand 

^^(^i'^)L,=.-.....-.„. =Y.^i{k,q)[kd,]^P\k,r^\ (52) 

1=0 

where q = k2 + ■ ■ ■ + km- The coefficients scale as hi oc (|g|/|/i;|)'. Therefore, 
since the Am grow as \k\'^^, we need to expand up to / < 2n to capture 
all terms that can potentially be affected by enhancement for large k. By 
collecting all terms arising from the /-th term in the Taylor expansion, and re- 
labelling the momenta, the ra-loop power spectrum for large k can be written 
in the form 

Pn-loop{k, V) ^ Y. I ^'^1 • • • ^'^- ^f"^^^' ^1' • • • ' ^-) 

X [kdkf P\k, v)PHki, V)--- P'-iK. V) , (53) 



^Thcsc quantities Pjj have nothing to do with the quantities Pat introduced above. It 
is always clear from the context to which quantity we refer. 
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up to terms that are suppressed by the hard external momentum, 0{ki/k). 
The Bi are hnear combinations of the Am for 1 < m < n -\- 1 multi- 
phed by the appropriate coefficients bi and we provide exphcit expressions 
in Appendix lAl Eq. (1821) . Therefore, without any cancellations, one would 
expect that Bi oc |A;p"~' for large k. The results from [11] show that the 
leading contribution cancels, which means that the enhancement is at most 
Bi oc |A;p'^"~^^"'. Our results of the previous section imply that this cancel- 
lation is also valid for the subleading terms and that 

Bl''\k, k,,..., kr,) ^ Ci''\k,, ...,kn) + 0{h/k) , (54) 

where B^ approaches a constant, denoted by C^ , for large k. 

6.2 Results 

An explicit analytical calculation showing the cancellation of k^- and sublead- 
ing /c^-contributions to the integrand kernels B^' up to two loops is presented 
in the Appendix El Explicit results for the asymptotic limit Eq. ( 15^ are given 
in Eq. fl68|) and Eq. ( 174]) for the one- and two-loop orders, respectively. In 
addition, we checked numerically for various momentum configurations that 
the cancellation of polynomially growing contributions and the asymptotic 
behavior Eq. flMjl is indeed correct up to four loop order, see Figs. [H] and O 
Using the analytical results Eq. fl^ and Eq. (1711) for the asymptotic 
values of the integrand kernels Bl it is possible to derive approximate ex- 
pressionqj (up to 0(l//c^) and subleading logarithmic corrections) for the 
power spectrum at large fc, 

Pi-ioop{k) ~ (1.14P^(A;) - 0.55A;9fcP^(A;) + Q.l[kdkfP\k)) af{k) , 
P2~ioopik) ~ (2.14P^(A;) - imkdkP'-ik) + 0.55[A;9fc]2p^(A;) 

- 0.082[kdkfP^{k) + 0m5[kdk]^P^{k))af{k) . (55) 

From these expressions, it is clear that the expansion parameter at large k is 
given by powers of (yf{k). We also see that even though all the functions B\ 
are of order unity, a logarithmic enhancement in the final power spectrum 
can arise depending on the precise shape of the initial power spectrum. For 



^We do not write explicitly the dependence on rj in the rest of this section, since it is 
trivial to retrieve. 
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Figure 8: Cancellation of n-loop contributions to the integrand kernels B\ of 
the power spectrum for / = 0. Shown are individual contributions Pjj (dashed 
lines) which grow like A:^" (indicated by the black dotted lines). Their sum yields 
the kernels -B; (red solid line), which approach a constant value as predicted in 



Eq. (j54p . For comparison, the grey lines show the analytical results (j68p . ()74p 
for the asymptotic value at one and two-loop. For the plot we have chosen ki = 
0.02, k2 = 0.03, k^ = 0.015, k^ = 0.019. For one-, two- and three-loop the kernels 
are averaged over the angles, see (I8ip . The four- loop case corresponds to (I82p . 



the spectrum (H^ one finds <jf{k) oc \v?{k) for large k. This is discussed in 
section [71 Additional enhancement can potentially come from the fact that 
the derivatives {kdk)^P^{k) are larger than the power spectrum P^{k) itself. 
Before moving to higher loops, it is interesting to remind the behavior 
at small k where the linear approximation is expected to work rather well. 
To understand the size of loop corrections parametrically, notice that in the 
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Figure 9: Asymptotic beliavior of n-loop integrand kernels B^ of the power 
spectrum for < ^ < 2n + 1 (here only the sum of all individual contributions to 
each kernel is shown). The kernels with / < 2n approach a constant value for large 
k as predicted in Eq. (|54p . and go to zero like 0{kf /k"^) for / > 2n. Parameters are 
chosen as in Fig. [H The grey lines show the analytical large- A; results ([68]) . (j7i|) 
for comparison. For / > we show only the range where none of the Heaviside 
functions contained in the integrand (I79p is zero. 



limit fc — ;■ 0, the dominating contributions to the n-loop corrections arise 
from the diagrams usually denoted by Pi,2n+i and scale as oc k'^P^ik) (see 
Appendix |B]). The behavior of this quantity is discussed in detail in [28] . 
We do not have much to add to this discussion but we want to identify the 
expansion parameter in this regime and compare it to the large k limit. For 
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the one- and two-loop expression^lj, we find 

fil Air f°° 61 

44764 Att 1"°° 

P2-ioo,{k) ^ "li^^'^ (^)y/ dqP\Q)A<l), (56) 

where we introduced the (dimensionless) function 

J{q) = An r dpp'g{p/q)P\p) . (57) 

An explicit form for g{x) is given in Appendix [Bl The relevant aspects are 
that it is a smooth function satisfying p^ g{p/q) = q^ g{q/p)- This means 
that g is almost constant as long as p and q are not of the same order and 
tends to zero for large argument. Thus, it further cuts-off the integrals of 
the form (l57|l at large value of the integration variable. We find (^(O) = 1.54, 
g{l) = 1 and g{oo) — t- 0. Thus, a good approximation for J(g) in the 
regime of integration of (1571) is given by taking g{p/q) ~ 1, which means 
that J{q) ~ crfiq)- Therefore, the relative importance of the two- loop results 
with respect to the one- loop case is again related to the scale ai{q). Notice, 
however, that in the low- A; case this function is integrated over (cf. fl57|l ). 

The previous calculations confirming the results of section \5\ can be ex- 
tended to higher loops. Even if the analytical calculations are very cumber- 
some in this case, one can still find general arguments about the different 
behaviors that may be checked with numerical computations. For the large 
k regime, as long as k is larger than the momenta one integrates over in 
(1^ . the results of section O imply a roughly constant Bi in (155]) . If one 
of the momenta ki becomes larger than k, additional polynomial suppression 
oc k'^/kf arises in the functions Bi what renders the integrals in internal 
momenta finite. Therefore, k effectively acts as a UV cutoff for the ki inte- 
grations, and one can estimate the integral to contain at n-loop contributions 
of order (/ < 2n) 

Pn-ioop 3 {[kdkfP'-ik)) af^ik) . (58) 

Subleading logarithms depend on how the integrals are precisely cut off and 
can give sizable corrections. 



^'^In the analytic calculations wc assume a flat matter-dominated cosmology and only 
consider the growing mode. 
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Figure 10: Linear power spectrum and the one- and two-loop corrections in stan- 
dard perturbation theory at z = 0. The hnear spectrum corresponds to a KCDM 
cosmology with WMAP5 parameters [26]. The blue lines show the asymptotic 
behavior at large k (see Eq. ([55]) ). and the red the one at small k (see Eq. (j56|) ). 



Concerning the low k case, from the symmetries of the integrand in 
and its behavior at large momenta, one expects the power spectrum to behave 



as 



n—loop\ 



(k) oc eP'^ik) / dqP'-iq) af"-2(g) . 



(59) 



As in the previous case, the important expansion parameter for perturbation 
theory is related to the quantity af. 

We have checked the previous asymptotic behavior by computing the 
power spectrum numerically up to two loops in an Einstein-de Sitter cosmol- 
ogy (taking into account only the growing mode). The results are displayed 
in Fig. [TOl We cross-checked our numerical results for the power spectrum 
with the RegPT code [26j for momenta where the latter is available. As 
can be seen in Fig. [TOl the asymptotic expressions f l56|) and f l55|) agree with 
the full one- and two-loop results rather accurately. The constant offset in 
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the two-loop case at very large k is expected because (jSSD captures only the 
leading logarithmic behavior. It is remarkable that even for very small mo- 
menta, the two-loop contribution is only mildly suppressed compared to the 
one-loop contribution. Naively, one might expect a suppression by a power 
/cqCt^ ~ (9(10~^), where ko is the position of the maximum of P^{k). How- 
ever, the integral over J{q) ~ o"f (g) in fISBl) is sensitive to the power spectrum 
at smaller scales, and yields only a mild suppression with respect to the one- 
loop correction. According to (159|) . the three- loop result should even exceed 
the one- loop contribution at asymptotically small momenta for z = 0. Ac- 
tually, we checked numerically for a few values of the momentum that this 
is indeed the case, and it is also in accordance with the findings of [2S]. The 
same is true for the large k regime. Here, the logarithmic dependence in 
crf{k) is clearly seen in the numerical results. 

7 Discussion 

Cosmological perturbation theory provides a very successful framework to 
understand the gravitational clustering in the Universe, responsible for its 
structure at the largest scales. Despite of this success, its range of validity (in 
the sense of convergence of the perturbative series to the non-linear solution) 
is limited due to the growth of the importance of the non-linear contributions 
with time. To devise methods that deal with these non-linear corrections it 
is crucial to understand which effects are behind the failure of perturbation 
theory. This is even more important when one realizes that a large amount 
of information about cosmological evolution lies at scales close to the linear 
scales where these methods may be very useful. An example is provided by 
the baryonic acoustic oscillations (BAO) at low redshift [3| |29| [30]. 

A first look at the structure of the standard formulation of cosmological 
perturbation theory (SPT) singles out the functions 

k'aj{A,v) = ^-y dqP'^iQ^v), (60) 

.A 

afiA^T]) = 4vry dqq'P'^iq^v), (61) 

as responsible for the failure of the linearized approach for scale k, once one 
of them becomes big. The first quantity, which is strongly k dependent, is 
small at the peak of the power spectrum ko, o'd^o — 0.0135 at 2 = and 
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ko ^ 0.02 Mpc~^ (see (HU]) ). This number grows very fast with momenta, 
which is a consequence of the enhancement coming from the vertex 7jjfc for a 
large hierarchy between the momenta of the two incoming modes, see (|8]). At 
n-loop order, this can potentially lead to corrections to the power spectrum 
that scale at large /c as oc {k'^a'^)"'P^{k,r]). It is well known that these 
leading soft corrections cancel when summing over all n-loop diagrams [TT] . 
However, there are also subleading soft corrections, growing like k"^ at two- 
loop, like k"^ and k'^ at three-loop, etc. In our analysis we showed that, 
when summing over all n-loop diagrams all polynomially growing corrections 
oc fc^™ with 1 < m < n cancel in the limit of large k. This cancellation 
had not been proven before to the best of our knowledge. Remarkably, the 
same cancellation happens for any hard skeleton corresponding to n-point 
correlation functions evaluated at equal-time. While the general proof relies 
on the eikonal approximation, we also checked explicitly that the subleading 
fc^-terms cancel at two-loop, and furthermore checked numerically that all 
subleading /^^'"-terms cancel up to four loop order. The cancellation of the 
leading soft corrections was related to Galilean invariance in [12] . Physically, 
it seems plausible that similar arguments could explain our results, though 
we are not aware of any explicit calculation. 

Our result also implies that in numerical calculations in SPT it is advan- 
tageous to sum over all relevant diagrams (and to symmetrize them appro- 
priately) before any integration is performed. In this way, the cancellation of 
different contributions occurs already on the level of the integrand and does 
not rely on the numerical accuracy of the integration. 

Even though the polynomially growing terms cancel, there remains a log- 
arithmic enhancement at large k for the matter power spectrum, see Eq. (153|) . 
We find that the leading logarithmic (LL) contributions at large k are given 
by0 



2n 



Pn-ioo,ik, v)^Yl '^^ [f^d^rP'ik, r^) a^-ik, rf) , (62) 

1=0 

with some coefficients q of order unity. This expression implies a growth 
with k of the expansion parameter af{k,rj) which is logarithmic. Note that 
af{k,'r]) also controls the loop expansion at the opposite limit of small k, see 
Eq. dSlD. 

The function af{k,ri) is not only sensitive to the high- momentum tail of 



^^ See Eq. ((55|) for the explicit expressions at one and two foops. 

27 



the spectrum (which makes it increase logarithmically) but for the realistic 
case f l49p it is also numerically rather large for small redshift 2; ~ 0, 

af{k,ri) c^ 0.15 (^—j (ln(e + 1.84/3 A;)) ^ (63) 

We find that this function is the true expansion parameter of standard per- 
turbation theory for the equal-time power spectrum. For /c ~ 1 Mpc~ it is 
crf{k,z = 0) ~ 3.42. This large value arises partially due to the logarithmic 
dependence but also because the initial power spectrum (149|) is parametri- 
cally enhanced by a factor (C/S^/cg)"^ ~ 16. According to our arguments a 
smaller value for ai would improve the convergence of linearized perturbation 
theory remarkably. We show this in Appendix O by using a fake initial power 
spectrum with the same ad as ACDM but smaller ai. The results are shown 
in Fig. [m which should be compared with Fig. [101 

Our conclusions about the matter power spectrum also hold for any cor- 
relation function at equal times. In these observables, the effects related to 
the scale a^ are absent, and the departure from the linear regime will be 
dominated by ai. The analysis is equally valid for arbitrary cosmologies, 
or even for departures from the standard equations, as the Zel'dovich ap- 
proximation. We would like to emphasize that this is not the case for other 
cosmological observables. For instance, the propagator as defined in ( 123|) 
is certainly affected by the enhancement of soft modes related to the scale 
(Trf. Its measurement in simulations confirms this behavior [1], and it would 
be very interesting to extract it from real data. Furthermore, correlations 
between different redshift bins are also used for lensing tomography [31]. 

All resummation schemes in the literature resum only certain subsets of 
diagrams of SPT. Notice that those subsets do not necessarily reproduce the 
cancellation we found for equal-time correlations. Thus, the non-linearity 
associated with the scale ad may be reintroduced as a purely spurious effect. 
We would like to emphasize that our results do not imply that standard per- 
turbation theory is superior to resummation schemes as e.g. RPT. It might 
well be that at intermediate scales, these schemes resum just the right sub- 
diagrams to lead to accurate results [32]. In addition, they are very useful 
to describe correlations at unequal times. Still, our analysis supports that 
at large momenta resummation schemes that involve only the scale ad can- 
not improve the determination of the equal-time correlators systematically. 
We hope that our results are helpful to identify approximation schemes that 
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respect the cancellation of soft corrections and, at least partially, resum cor- 
rections related to the scale ai. 
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A The power spectrum in the hard regime 



In this appendixo we discuss the asymptotic behavior of loop corrections 
to the power spectrum at large wavenumbers k (small scales). It is possible 
to rewrite the loop contributions such that the cancellation of polynomially 
growing corrections oc k"'P^{k), with n > 0, is manifest. We first demon- 
strate this for the one- loop corrections Pi^-ioop = 2Pi3 + P22, with 

Pi3(A;) = SP'^ik) fE!{k,q,-q)P\q), (64) 

Jq 

P22{k) = 2 /"k(g,^-g-)]'p^(g)P^(|^-gl), 

where P^(gi, . . . , g^) are the symmetrized PT kernels entering at the dif- 
ferent orders of SPT (see e.g. [2]), / = J d^q and q = \q\. For large k 
the two contributions asymptotically grow as P13 -^ —k'^a'^P^{k)/2 and 
P22 — )■ k'^cr'^/2P^{k), and it is evident that the quadratically growing cor- 
rections cancel [33] • At two-loop, P2-100P = 2P15 + 2P24 + P33, the individual 



^^To improve readability we will not write explicitly the dependence on 77 in the ap- 
pendices, since it is trivial to retrieve. In addition we display arrows on three vectors to 
faciliatc discrimination from absolute values. 
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contributions 



P,5{k) = 15P\k) f F^{k,p,-p,q,-q)P\p)P\q), (65) 

P,,{k) = 12 f F^{q,k-q)Fi{q,k-q,p,-p)P\p)P'^{q)PWk-(^), 
Jp,q 

P,,{k) = 9 f F^{k,p,-p}F.^{k,q,-q)P''{p)P'^{q) 

+ 6 / Fi{p,q,k-p-q)^P\p)P'^{q)PWk-p-q\), 
Jp,q 

grow asymptotically as k^a'^P^{k). Again, it is straightforward to check that 
the k"^ terms cancel each other. However, there exist also subleading terms 
that grow as fc^. In order to show that they also cancel in the sum of all 
two-loop contributions, it is convenient to rearrange the various terms. To 
demonstrate this, we first discuss an analogous rearrangement for the one- 
loop terms, 

P,L{k) = P'^ik) f B'i\k,q)P\q) + P^^{k), (66) 

Jq 

where {dVlq = sin OqdOqdipq) 

Bi'\k,q) ^ ^ J dnq{QF^{k,q,-q)+AF^{q,k-q)') , (67) 

P22{k) ^ AJF^{q,k-q)^P\q)(Q{\k-q\-q)PWk-q\)-P\k)). 

In the latter contribution we subtracted a term oc P^{k), such that the 

difference in the bracket scales as (|f -|- 0{q^ / k'^))kdkP^ {k) for large k. The 
linear term oc q/k vanishes when integrating, such that the bracket leads 
to a g^/A;^ suppression that compensates the quadratic enhancement coming 
from the kernel (-^2'*)^ ~ {k(i/(fY /^- Note that we also inserted a Heaviside 
function compared to P22 and multiplied by two, which does not change P22 
due to the symmetry g — )■ k — q. Using the explicit expressions for the kernels, 
it is also easy to check that quadratically growing terms ~ k"^ /q^ cancel in the 
combination B^ after averaging over angles (or alternatively symmetrizing 
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the integrand w.r.t. q ^ — g). The exphcit result can be easily calculated, 



B^^\k,q) 



625 llk'^ 



523k^ 



q^ 



882 29V 1764g2 I2k^ 



k^-q' 
q2 



llk^ 



9k 8q 

+ 



1176^3 112g A9k A8k^ 



r \i^(^~^ 



k + q 



2205 ^ ^"' ' ' 



for k ^ oo . 



(68) 



A similar rearrangement can be done for the two-loop contributions to the 
power spectrum. In order to extract subleading fc^-terms apart from the 
leading /c^, we have to re-shuffle terms proportional to the first and second 
derivatives of the power spectrum (' = d/d\nk), 



P^{\k - q\) = P'^ik) + bi{k, qjP"^ (k) + b2ik, q)P'^ (k) + . 



(69) 



where 



bi{k,q) 



In 



\k-^ 



(70) 



After this rearrangement, the two- loop contribution can be rewritten as 
P2Lik) = P\k) I B^;}\k,p,q)P\q)P\p) 

+ P'^'{k) f i?f)(A:,p,g)P^(g)P^(p) 



?(2), 



+ P^ (A;) / B','>{k,p,q)P'^{q)P^{p) + 2P2,{k) 
Jp,i 

+ P33{k), 



(71) 



where 



?(2) 



Br>{k,p,q) 



J d^g j dnp(30F^ik,p, -p, q, -q)5, 

+ 9F^{k,p,-p)F^{k,q,-q)6io 

+ 18F3^(p, q,k-p- q)\{k,p + q) 

+ 24F2'(g, k - q)Fl{q, k - g, p, -p)hi{k, q) 

+ 24F2^(p, k - p)Fl{p, k-p, q, -q)kik,p) 



10 



(72) 
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for Z = 0, 1, 2 with bo = 1, and 

P2,{k) = 24 f F^{q,k-q)F!{q,k-q,p,-p)P'^{p)P\q)(Q{\k-(^-q) 
Jp,q 

xPWk-c^)- P\k) - h{k, qjP'^'ik) - hik, q)P'''\k)) , 

Pssik) = 18 f F^{p,q,k-p-q)'P'^{p)P''{q)(e{\k-p-q]-q) 
Jp,q 

X Q{\k - p - q\ - p)P^{\k - p - q\) - P^{k) - h{k,p + q)P^' {k) 
-b2ik,p + q}P'^"ik)). (73) 



The subtracted terms inside the brackets in the last two expressions are 
constructed such that, after angular integration, they compensate the terms 
growing as k^ and fc^ coming from the PT kernels for /c — )■ oo. Thus, it 
remains to be shown that the functions B^ in ( ITT]) do not grow with k. For 
that purpose, we expand the integrand for large k and use that the angular 
integration can equivalently be performed over k and p while keeping the 
direction of q fixed along the 2-axis (see e.g. [I2|)- We find that indeed the 
k"^ and k"^ terms cancel. For the leading contributions in the limit /c — )■ oo 
(which are oc /c°) we find 

„f2)., ^ 271133 /,, , , 595172936\ ,^,,,,2^ 

B\:Hk,p,q) -^ fip q) + +0(1 k^), 

V ,y,Hj 92702610V 813399 J ^' '' 

^(2),, , 5377 / , ,, 308888266\ ,^, ,,2, 

«'(*•"• ') - 12677280 V''^'"^ - ^065^) + ^'^/* * ' 

Bf («:.?. 9) ^ ^ {f(P/l) - ^) + 0(1/*') , 
where 

The results for / = 3, 4 are shown for later use, and Bj^ = 0{\/k'^) for / > 5. 
For completeness, we also quote the large-/;; limits for the higher-derivative 
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contributions at one-loop which are given by 

Bi'\k,q)^-'^^, Bi'\k,q)^^, Bl'^,{k,q) ^ 0(1/^). (76) 

The rearrangement presented above is useful to show that polynomially 
growing corrections cancel out. The leading correction for large k is thus 
a logarithmic one. It is possible to extract an analytic expression for the 
leading logarithmic corrections by performing a similar rearrangement, but 
including terms up to [kdh^P^ at one- loop and up to the fourth derivative 
[kdkl'^P^ ^^ two-loop. This yields approximate expressions for large k 

Pi-ioor>{k) -^ (|^^"(^) - '%kd,P\k) + ^^[kd,fP\k)^ af{k) , (77) 

4 pk pk 

P2-ioop{k) ^ (AnfJ^lf'dkfP^k) / dqq^P\q) / dpp'P\p)B^;'\k,p,q) 



1=0 "^0 "^0 

/l490372537„^,,, 7719787, ^ „i,, , 112327„ ^ .3 ^^ ,, , 

~ P (k) kdkP ik) H \kdk\ P ik) 

y 695269575 ^ ' 4753980 ^ ' 205800 ^ ^ ^ ' 

1207 1 \ 

~ 14700 f^^'^]'^'^^^ + ^[^^fc]'^"^!^) j^f(^) • (78) 

Here we assumed that the main contribution to the integrals in (!78|) comes 
from the integration range p,q < k, which given the form (H9|) is valid at 
high k up to logarithmic corrections. Besides, in the last line we substituted 
f{p/q) ~ 1, which is valid for p ^ g (or, equivalently, q <^p). One can easily 
get convinced that this estimate is adequate for /(x) given by fl75|) . 

It is straightforward to extend the rearrangement to n loops. The n-loop 
contribution given in Eq. f lSTjl can be rewritten by integrating over ki. Due 
to the symmetry in the momenta ki, . . . km, one may restrict the range of 
integration to the case where \ki\ is larger than any of the \ki\ for 1 < i < m, 
and multiply by a factor m, 



n~loop 



n+1 n 

{k,ri) = ^m j d^k2...d^kn+iA^Z\ki.---.kn+i) 

m=l "^ 



y<Qi\kl\ - \k2\) X ■ ■ ■ X Q{\k,\ - \kra\) 

I ki=k—k2 k„ 



xP\k,,v) . ..P\k^^,,rj)\,,r, . (79) 
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The integration kernels in the notation of Eq. (151 p are given by 

n—m+l 



yl£'(fc,...,*„+i) 



ni=0 



{2nL + m)\{2nji + m)\ 



^ -^2ni+m(^l5 






^ -^2nfl+m(^l5 • • • , «^' 



I'^miPli Pli ■ ■ ■ iPnx^i Pn_ty 



(80) 



where n^ 



n 



m 



riL, Pi 



k 



m+i 1 Hi 



k 



m+riL+i- 



In the usual Pjj 



notation, a given term in the sum in (|80|1 contributes to -P2ni,+m,2n«+m- Dia- 
grammatically, a given summand corresponds to a diagram with two 'blobs' 
that are connected by m lines, and have ul and rtn lines connected to them- 
selves. By adding and subtracting the terms obtained from a Taylor expan- 
sion of P^iki^Tjo) around k up to order 2n, one obtains the rearrangement 
analogous to the one- and two-loop case discussed above. The coefficients of 
the terms containing the Taylor-expanded power spectrum in ( !53|) after the 
rearrangement are 



Bt'\kM.. 

where 

Bt'\kM, 



) i^n) 



(47r)" 



dVlk^ 



An) 



drtk„Bl '{k,ki,...,kr, 



(81) 



; kn) 



AP{k,ki,...,kn)Sio 

n+1 
m=2 



km-1^ k 



1) 



; kn) 



xbi{k,ki + 



k 



m — ly 



symm 



(82) 



Here, the right-hand side is to be fully symmetrized w.r.t. permuting the ki 
and w.r.t. inverting the sign ki — ?■ —ki of each momentum. We note that in 
our numerical results (up to four loop order) we find that the cancellation of 
terms growing with k can even be observed at the level of the non-averaged 
expression. 
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B The power spectrum in the soft regime 

In the limit /c — > 0, the dominant contribution to the n-loop correction to 
the power spectrum is given by [211 ^^ 



n~loop 



(fc) ^ 2 (2n + 1)!! P^(A;) /"■■■/" Fi^^,{k, gl, -gl, . . . , q^, -Qn) 

•J qi J Qn 

xP^(gi)---P^(gn). (83) 



Due to the well-known property F2nj^i oc k'^ of the PT kernels [2Z], these 
corrections scale as k'^P^{k). All other loop corrections would lead to terms 
scaling as /c^ or k'^P^{k), respectively, that are parametrically smaller for 
fc — )■ for a power spectrum P^{k) oc /c"" with n^ < 1. The function 
appearing in the two-loop correction ( 15B|1 is given by 



g(x) = — (x^ + 1) (I28258a;^-5760(a;^ + l)-13605(x^ + x^)) 

I79056x° ' 



— ix^ - If (384(x^ + 1) + 2699x2) In (- — ^ 
Ax \x + \ 



^4) 



C Convergence and fictitious power spectra 

In this section, we present a two-loop result for a fictitious power spectrum 
with better convergence behavior than fj^9|) . Let us consider an initial power 
spectrum of the form 

Po{k) oc ^^ . (85) 

This spectrum is chosen such that the expansion parameter of standard per- 
turbation theory af is not sensitive to the high momentum part of the spec- 
trum. Hence af is small as long as k'^a'^ is much smaller than unity. Fig. [TT] 
shows a power spectrum with k^a^ ~ 0.015 and af ^ 0.02. Even though /cqO"^ 
is as large as in the physical case at 2 ~ 0, the two-loop result is well below 
the one-loop contribution and standard perturbation theory is expected to 
converge. Compare with Fig. [TDl 
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Figure 11: Example for a power spectrum where af is not sensitive to the 
high momentum region. The different hues are the hnear, one-loop and two- loop 
contributions. The corresponding parameters are /cqct^ ~ 0.015 and af ~ 0.02. 
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